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ABSTRACT 

We present the wide-field imaging and polarimetry at v = 20 GHz of seven of the 
most extended, bright (Stotal 0.50 Jy), high-frequency selected radio sources in the 
southern sky with declinations S < —30°. Accompanying the data are brief reviews of 
the literature for each source. The results presented here aid in the statistical com- 
pleteness of the Aust ralia Telescope 20 GHz Survey (AT20G) Bright Source Sample 
l|Massardi et al.ll2008t Burke-Spolaor et al., in prep). The data is of crucial interest for 
future CMB missions as a collection of information about candidate calibrator sources. 
We were able to obtain data for seven of the nine sources identified by our selection 
criteria. We report that Pictor A is thus far the best extragalactic calibrator candidate 
for the Low Frequency Instrument (LFI) of the Planck ESA mission due to its high 
level of integrated polarized flux density (~ 0.50 ±0.06 Jy) on a scale of 10 arcminutes. 
Six out of the seven sources have a clearly detected compact radio core in our images, 
with either a null detection or less than 2 percent detection of polarised emission from 
the nuclei. Most sources with detected jets have magnetic field alignments running in 
a longitudinal configuration, however PKS 1333-33 exhibits transverse fields and an 
orthogonal change in field geometry from nucleus to jets. 

Key words: polarization - radio continuum: galaxies - galaxies: magnetic fields 
- galaxies: active - cosmology: cosmic microwave background - galaxies: individual: 
Centaurus A, Centaurus B, Pictor A, PKS 0131-36, PKS 1333-33, PKS 1610-60, PKS 
2153-69, PKS 2356-61 



1 INTRODUCTION 

The Australia Telescope 20 GHz survey ( AT20G) is a, blin d 
20 GHz survey of the whole southern sky (jRicci et al.ll20o4 ) . 
The project outlined in this paper complements two sub- 
projects of the AT20G Survey; the fir st is the AT20G Bright 
Source Sample (Massard i et al.ll2008l . BSS), and the second 
a high-sensitivity polarisation study of a subset of the BSS 
(Burke et al., in prep). The collective goal of these sister 
projects was to analyse the total intensity and polarisation 
statistical properties for a flux density limited sample of 
sources selected at 20 GHz. The BSS used sources south of 
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5 — —15°, while the Burke sample focussed on sources south 
of S = —30°. Both samples included all extragalactic AT20G 
objects which were brighter than 0.50 Jy, and are primarily 
comprised of a population of compact objects, with 85 per- 
cent showing a f lat spectrum (a > — 0.5; S v oc v a ) at high 
radio frequency l|Massardi et al.ll2008h . 

This paper enhances the two projects by providing flux 
density and polarisation data for objects of angular extent 
beyond approximately 2.4 arcminutes, for which accurate 
data could not be obtained using a single interferometric 
synthesis field. This paper includes sources only sources 
south of S = -30°. 

Information about nearby, bright sources is of im- 
portance to observations of the Cosmic Microwave Back- 
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ground (CMB), in particular the current focus on detect- 
ing anisotropics in the polarised CMB signature. Thomp- 
son scattering on the Last Scattering Surface produces lin- 
ear polarisation and curl-free patterns termed "E-modes," 
which on small angular scales provide information about 
structures in the early Universe. Observations of compact 
CMB foreground sources are necessary to distinguish be- 
tween real E-modes and signal scattered by gravitational 
lensing of the CMB on point source angular scales. Addi- 
tionally, the diffuse plasmas in bright, extended extragalac- 
tic sources c an form patterns which are termed "B-mode" 
polarisation (|Hu. Hedman. and Za ldarriaga 2003); B-mode 
signal in the radio lobes of extragalactic sources may be un- 
surprising and can be caused by complex magnetic fields 
in the sources and irregular substructures of different Fara- 
day rotation. However, the positive detection of a similar 
B-mode polarisation signal in the CMB is a major goal 
of increasingly sensitive measurements of CMB anisotropy. 
A B-mode signal in the CMB is predicted to occur only 
by CMB models with tensor (gravitational wave) fluctua- 
tions; these are in turn tied to the rapid period of expansion 
that may have occured prior to the recombination epoch. 
The detection of B-mode CMB anisotropics will give insight 
into the energy scales that occurred directly prior to infla- 
tion and would suggest the presence of gravitational waves. 
However, the amplitude of the cosmological B-mode sig- 
nals is extremely small, and will require careful attention 
to foreground signals that might mimic it. On small angu- 
lar scales where extended and complex foreground sources 
(such as radio galaxies) are resolved, removing contamina- 
tion due to extragalactic sources becomes particularly im- 
portant. For faint, confused, and unresolved extragalactic 
sources in CMB data, a statistical correction may become 
necessary to calibrate and interpret the results of a cosmic 
background B-mode analysis. 

Bright, highly polarised sources are also necessary 
for upcoming CMB instruments, in particular, the ESA 
Planck satellite mission's low frequency instrument (LFI), 
which will observe in bands between 30 and 70 GHz 
l|Bersanelli and Mandolcsi 200(J). Among the source types 
that have been pursued as potential flux density and po- 
larisation calibrators for the instrument are extragalactic 
sources; however, as the spectral behaviour of the bright ex- 
tragalactic population at high frequency is difficult to model, 
sources selected at lower frequency are not certain to be 
suitable calibrator can didates. With our obse rvations, along 
with the catalogues of iMassardi et ah I (|2008h and Burke et 
al. (in prep), we provide a list and homogeneous observations 
of bright and strongly polarised sources that are selected at 
a frequency closer to CMB observing frequencies than was 
previously available. 

In ij2] we outline the criteria used in the selection of 
our sources and summarise the literature accompanying each 
source. [JJ] details the data collection, calibration, and wide- 
field imaging process used to create our results. Finally, the 
data (full intensity and polarisation vector images and val- 
ues) are presented in Sjjjwith a discussion about the images, 
high frequency polarimetry, and the implications of our re- 
sults. [J5] provides a brief summary of our findings. 



2 SOURCE SELECTION AND LITERATURE 
REVIEW 

The AT20G team's confirmation followup observations to 
their initial transit scan survey provided a list of sources 
detected at 20 GHz, from which Burke et al. (in prep) se- 
lected a sample of sources at declinations —30° > S whose 
total intensity exceeded 0.50 Jy. Pr eliminary res u lts of an 
AT20G pilot survey are available in ISadler et al.l (|2006l ): a 
description of the blind scan and f ollowup observa t ions o f 
the AT20G survey are available in IMassardi et all (2008). 
The 2.4 arcminute primary beam FWHM and the 80 to 
200 metre interferometer spacings used in the AT20G fol- 
lowup did not provide accurate flux density measurements 
for the most extended sources, and furthermore left cases 
in which sources were either: 1) Fully resolved (and there- 
fore undetected) by the 60 metre shortest antenna spac- 
ing of the AT20G followup, or 2) Had compact compo- 
nents (hotspots, cores) which were detected as separate, 
<0.50 Jy sources in the followup. T he Sydney Univer - 
sity Molonglo Sky Survey (SUMSS) dMauch et al.1 I2003T) 
and the Parkes-MIT-NRAO (PMN) South ern (|Wright et"al] 
11994 [Condon. Griffith, and Wrightl 1 1993ft and Zenith Sur- 
veys l|Wright et al.l Il996h have a higher sensitivity to ex- 
tended radio sources. This is due to their more complete 
Fourier coverage for short baselines, and due to the lower 
observing frequency of the two surveys which allowed higher 
sensitivity to diffuse, steeper spectrum radio components 
(SUMSS at 0.843, PMN at 5 GHz). 

To recover highly resolved 20 GHz sources, we cross- 
referenced multiple-component AT20G sources against 
SUMSS images. This indicated AT20G sources which ap- 
peared to be compact components of one larger, extended 
object. We then ran a search for sources which are highly 
extended and thus may have been completely resolved out 
by the AT20G obervations. We expect any extended source 
with S 2 o > 0.5 Jy to have a 5 GHz flux of S 5 > 0.9 Jy if it 
has a spectral index of less than -0.5, as would be expected 
by a conservative estimate for the spectral index of diffuse 
emission. We thus compiled a list of the sources in the PMN 
catalogue which are flagged as extended and have an inte- 
grated 5 GHz intensity greater than 900 mjy. By including 
these objects and those selected at 20 GHz, this sample has 
a nominal flux limit for extended sources with S20 > 0.50 Jy. 
However, due to the spatial resolution limit of PMN and our 
decision to include only sources flagged as extended in PMN, 
the sample is not complete for sources of scales between 2.4 
and 4.2 arcminutes. The number of missing sources in this 
range is estimated to be between 4 and 10 based on source 
count extrapolations. 

Sources which were associated with the Large and Small 
Magellanic Clouds (LMC, SMC) or in the Galactic plane 
(|6| < 1.5) were excluded from the sample, as were supernova 
remnants with \b\ > 1.5. 

In summation, a source was included if it: 

• is in the region -30° < 5 < -70° 

• is extragalactic (not including LMC, SMC). 

• was detected as a multiple-component (hotspots/core) 
source within the AT20G followup beam with a total flux of 
Stat > 0.50Jy. 

or 

• is marked in PMN as extended, with 5*5 > 0.90 Jy. 
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Table 1. The nine sources identified by our selection criteria. Columns 1 and 2 give the common source name and the WMAP ID, 
respectively. Columns 3, 4 give the radio centroid position. Column 5 gives the source redshift; for redshift references, see text. Columns 
6, 7 Give the angular size and calculated linear size based on a flat, = 0.7, f2 m = 0.3, Ho = 72 km/s/Mpc cosmology. Column 8 gives 
relevant details of the observations (including number of telescope pointings for observed sources). 



If a source obeyed all the above criteria, it was included in 
our sample. 

The resulting list from the searching process consisted 
of nine objects. However, we did not observe Fornax A due 
to its highly diffuse emission, most of which would remain 
undetected even in ATCA's most compact configuration 
(furthermore, attempts to detect the nucleus of the galaxy 
at 20 GHz during unscheduled ATCA telescope time 
were unsuccessful down to a few millijansky) . For some 
sources, only subregions containing compact structure were 
observed. The full nine sources detected by our criteria are 
listed in table [T] which provides information on what part 
of the source was observed (if applicable). 

Below, the radio morphology of each source is described, 
and significant literature and interesting properties for each 
object are summarised. 

PKS 0131-36 (NGC612), figure H An asym- 
metric, Faranoff-Reilly II source with a bright kn ot in 
one jet at low radio frequency (|Mauch et al l 12003). As- 
sociated at optic al wavelengths with a magnitude 1 5 
SO-type galaxy (iBolton. Gardener, and Mackevl 1 1964 ). 
IWesterlund and Smith! l|l96rj ) note a radio-optical mor- 
phology similar to Centaurus A, with a dust lane crossing 
the core and an asymmetric positioning of the core along 
the axis of the two radio lobes. Recently, a study of HI 
gas in the galaxy demonstrated the existence of a huge, 
140 kpc disk along NGC612's dust lane and a faint bridge 
of HI emission to the nearby galaxy NGC619, indicating 
a possibly recent o r ongoing interaction in this system 
l|Emonts et al.ll2008l) . The source re sides near the centre of 
a cluster of galaxi es (|Schilizzilll975l) . Redshift measured by 
lEkers et all (|l978l ). 



Fornax A (NGC1316): Fornax A is a highly diffuse, 



symmetrical double-lobed galaxy extending approximately 
one degree along its major axis. So urce is assoc i ated with 
the p eculiar SO galaxy NGC1316 l|Millsl Il954l ; iFergusonl 
Il989l ) . The source resides in the Fornax cluster. The high 
resolution observations of iGeldzahler and Fomalontl (|l984l ) 
across 1.5, 4.9, and 15.0 GHz indicate a steep-spectrum 
core and nonlinear jets. T he large-scal e stru cture of the 
radio source was studied by lEkers et alj (|1983l ). who noted 
a bridge of emission extending between the lobes that 
is displaced from NGC1316's centre; this morphological 
feature and the complicated dynamics of the galaxy lead to 
their description o f a formation scenario involving multiple 
galaxy collisions. ICaon. Capaccioli. and D'Onofricl (|l994T ) 
further suggest Fornax A's host as a merger remnant, noting 
tidal tails, filaments, dust lanes, and a low velocity disper- 
sion. Although too diffuse to make useful measurements 
with our telescope configuration, the source was detected 
in WMAP at 23 G Hz with a flux of S23 = 11.7 ± 0.1 
l|Bennett et al.ll2003l). The redshift given in the table is as 
measured bv lLonghetti et al.l l|l998h . 

Pictor A, figure [3] A classic double radio galaxy 
with si gnificant hotsp ots, identified both as a D type 
galaxy (|Schilizzil Il975l ) and as an N typ e galax y with 
strong X-ray emission by Mar shall et al.l (|l978l ). It is 
known to be hosted by a smooth, e ~ 0.7 elliptical 
galaxy wit h a sharp nucleus and little to no fil amentary 
structure (Danzige r. Fosburv. and Penstonl 1 19771) . Exten- 
sive observations and analysis have been done of this 
source. These include studie s of it s rot ation measur e, for 
instance iBerge and Seielstadl l|l967h and iHavesl (fl975) give 
RM = 46 ± 1, 6> = 104 ± 4 and RM = 51 ± 4, 6> = 100 ± 6, 
respectively, where 9o is the calculated zero-wavelength 
polarization angle. In addition, there have been high- 
resolution observations (|Preuss and Fosburvl i|l983h note 
a radio nucleus of parsec scale) and central black hole 
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mass estimations (|Lewis and Eracleousl l|200d ) estimate 
Mbh = 4 ±2[10 7 Mq] using stellar velocity dispersion) .The 
western radio lobe contains a bright hotspot that shows 
highly polarised (optical polarisation degree ~55 percent) 
optical synch rotron emission and coinc id ent X -ray signal 
from which iRoeser and Mei scnhcimcr (1987) conclude 
pure synchrotron emission from radio to X-ray wave- 
leng ths. A deep optical study o f the w estern hotspot 
by iThomson. Crane, and Mackavl |l995) confirm this 
conclusion. The study also notes strong polarisation of 
the hotspot with a magnetic field oriented perpendicu- 
larly to the jet axi s. A thorough X-ray study o f Picto r 
A was done by IWilson. Young, and Shopbelll i|200ll ). 
noting significant emission in the core, bright western 
hotspot, and along a jet extending from the core to 
the hotspot. High quality imaging of Pictor A is avai l- 



able at vari o us wa velength s (e.g. | Tingav et al 



Simkin et al.l (Tl999T); iMalkan, Goriian. and Tarn 



2000) 



1998 



Perlev. Roeser. and Mcis enheimerl (|1997T )). Redshift mea- 
sured by ISchmidtl (|l965T )~ 



Centaurus A (NGC5128), figure [4] A very complex 
and highly extended radio source associated with a spherical 
galaxy crossed by large, absorptive dust lanes. Being the 
nearest and one of the first radio galaxies identified, the 
literature is incredibly extensive and too vast to summarise 
here. A seminal study of the large scale radio structure 
of the source from 500 MHz to 5 GHz was made by 
ICooper. Price, and Colel (ll965T): the optical identification 
was made bv lBolton. Stanley, and Sled (|l949Tl . The redshift 



measurement is available from Graham! (1973) 



PKS 1333-33 (IC4296), figure O A "triple" 
radio source called thus for its bright inner jets and diffuse 
outer lobes that are in alignment along the jet axes. The 
major lobes extend approximately 30 arcminutes and 
the inner lobes a r e within 2 arcminutes of the centroid 
(|Goss et all Il977l ; Ijones and McAdaml Il992f) . The radio 
source is associated with the El galaxy IC 4296, which is a 
companion to IC 429 9 and is centrally located in the galax; 



alaxy 
1985; 



cluster Abell 3565 ijYounis. Meaburn. and Stewart] _ 
iHuchra and Gellerll 19821). The host galaxy is claimed to have 
a du sty disk by Colbert. Mulchaev. and Zabludofi ((2001 



and ISchmitt et afl~ ( 2002) : however, Sadler and Gerhard 
(|l985h and iMichardl (ll998Tl saw no sig nifica nt dust 



features. Killeen. Bicknell. and Ekersl (|l986T ) and 



iKilleen. Bicknell. and Carter! ( 19861 ) provide thorough 



multi-wavelength studies of the radio polarisation and 
rotation measures, and studies of the optical, infrared, 
and X-ray f orm of the source w i th an extensive analysis 
presented in IKilleen and Bicknelll (1 198ST ) . The redshift was 
most recently measured by Smith et alj (|2000l ). 



a spectral index of -0.73. Extrapolating their published 
measurements, we predict an integrated 20 GHz flux value 
of S20 ~ 1 4.6 Jy. The redshift quote d in the table was 
measured by I West and Tarenghil l|l989l ). 

PKS 1610-60, figure [7] The radio source is a 
wide angle head-tail galax y with bright inner lo bes and 
diffuse tails (mapp ed by jSchilizzi and McAdaml (1 19751) ; 
Jones and McAdaml (|l992T l: IChristiansen et al.l (1 19771 ); 
Gregory et al l d 19941) at 408. 843, 1415, and 4850 MHz, 



respectively). [Ekersl 1 19701) identified the source with an E3 
galaxy at J2000 a =16:15:1.7, S = -60:55:13.7. The galaxy 
is thought (but not yet proved) to b e a central member of 
the lo c al massive cluster Abell 3 6 27 llRobertson and Roachl 
Il99d ; H rown and Burns! Il99ll ; iBoehringer et all ll996f T 
Another head-tail galaxy lies within one degree of PKS 
1610-60; I Jones and McAdaml (|l996l ) mapped the two 
sources at 1360 MHz, explaining PKS 1610-60's morphol- 
ogy (sharp kinks in the jets and the wide tail separation) 
through probable interaction with an intcrcluster medium 
and influence from subclu ster mergers. Most recently, a 
redshift was determined bv lWoudt et al.1 (|2004l ). 

PKS 2153-69, figure [5] Less than 3 arcminutes in 
extent, the radio structure of this Seyfert II galaxy consists 
of a core-lobe morphology. A presumably unre lated source 
lies approximately 2.5 a rcminutes to the east l|Ekerslll969l ; 
IChristiansen et aL 19771). The host galaxy is an early type o f 



class E0 (ISchilizzilll975l ; iMarenbach and Appenzellerlll982r i. 
and has been shown to have an extra-nu clear ionised gas 
cloud with which the jets are interacting (jTadhunter et all 
1 19881 ; iTingav et aT] ll996T ). Furt her deep radio and X -ray 
studies have been performed by iFosburv et alj (|l998l ) and 
lYoung et all (|2005l ). The galaxy was modell ed to have a 



long-p eriod (t ~ 1.8 x 10 6 y) jet precession by Lu and Zhou 
|2005l ). Redshift given by |Pa Costa et al.l (|l99ll ). 

PKS 2356-61, figure H A Faranoff-Riley II galaxy 
marked by four bright regions of emission that are 
slightly asymmetric about a core. In addition, there is 
a diffuse arm of emission extending south-west from 
the weste rn inner lobe (see e. g. low radio frequency im- 



Christiansen et al. (1 19771); Jones and McAdam 



Subr ahmanva n~Saripalli. and Hunsteadl (| 19961 ) 

iBurgess and Hunsteadl (I2006T I). Id entified with an E3 
galaxy by IWhiteoakl l| 19721) and ISchilizzil (|l975l ). Sus- 
pected to be a memb e r of g alaxy cluster SC 2357-61 
dWesterlund and Smith! Il966l ; iTeague. Carter, and Gray] 
199d). Source was detected as a hard X-ray emitter by 



Landi, Malizia, and Bassani (2005). Redshift measured by 



Lovedav et all i|l996l ) 



Centaurus B, figure [S] This is a double-lobed 
source spanning approximat ely 15 arcminutes, identified 
with a large, red E0 galaxy bv lLaustsen. Schuster, and Westl 
l|l977l ) . The galaxy is close to the galactic plane (b ~ 1.73) , 
and as noted by Ijones. Llovd. and McAdaml |2001), 
literature flux estimat es are often confused by g alactic 
foregrounds. However, Ijones. Llovd. and McAdaml (|200ll ) 
present a comprehensive multi-frequency study of the radio 
galaxy, finding a close power-law fit to the spectrum with 



3 OBSERVATIONS AND DATA REDUCTION 

All sources were observed with the Australia Telescope Com- 
pact Array (ATCA). The compact configuration H75 has 
five 22m dishes extending on perpendicular north-south and 



1 An unpublished image of this source is available from the ATNF 
at http:/ /www. narrabri.atnf.csiro.au/public/images/pks2356-61 
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east-west arms with a maximum baseline of approximately 
80 metres. The array was well-suited to obtain sufficient 
(u,v) coverage to sample and image large, diffuse structures 
at high radio frequencies over a short period of time. Each 
source was observed simultaneously at 16.7 and 19.4 GHz 
in two polarisations using linear feeds with 128 MHz band- 
width at each frequency. Observations took place on October 
1, 2006. 

The dish size and maximum baseline correspond to a 
field of view of 2.6 arcminutes and a resolution of 43 arcsec- 
onds at our mean observing frequency of 18.0 GHz. Because 
all the sources extended beyond the primary beam size of the 
telescope, each source was observed using ATCA's "mosaic 
mode" ; the mosaicing technique is commonly used for wide- 
field imaging at the ATCA. Overlapping, adjacent telescope 
pointings that were later jointly deconvolved as one obser- 
vation were taken in a hexagonal pattern with field centres 
spaced at approximately 1.5 arcminutes. The field spacings, 
given by the Nyquist sampling limit for a hexagonal grid, 



(where 6 is the angular size of pointing separation and 
D is the primary beam FWHM) ensure a constant noise 
level across overlapping regions of observation. Additionally, 
the joint processing of the hexagonally interleaved telescope 
pointings give the interferometer a sensitivity to large scale 
emission that is comparable to single-dish (zero-spacing) 
scales. Each mosaic subfield was observed for 40 seconds, 
and a secondary calibrator was observed before and after 
each mosaiced source. Calibrator and mosaic sets were per- 
formed four times for each source at a range of hour angles. 

All data were calibrated using the Miriad software pack- 
age (|Sault. Teuben. and Wrightlll993T ). Throughout calibra- 
tion the two IF bands were treated independently. After 
opacity corrections and xy feed phase difference solutions 
were applied, the bandpass was calibrated using a standard 
bandpass calibrator (PKS B1921-293), and a primary flux 
density scaling was solved for using a standard ATCA pri- 
mary flux calibrator, PKS B1934-638. These calibrators were 
both observed once only, directly prior to the observing run. 
Each secondary calibrator was used to simultaneously solve 
for antenna gain matrices, the residual xy-phase difference, 
and polarisation leakage terms (which result from slight im- 
perfections in feed alignment). Solutions from this process 
were then applied to each subfield. In the case of Centaurus 
A, additional self-calibration with a single Gaussian compo- 
nent model was needed for fields containing the bright core 
of the source. 

Due to current limitations and software incompatibility 
in some of the deconvolution algorithms, our measurement 
of total polarised flux and the images published in this work 
show the results of Miriad's joint deconvolution polarisa- 
tion maximum entropy method algorithm (pmosmem), while 
the total intensity measurements are extracted from a new 
multi-scale clean deconvolution algorithm implemented in 
the Common Ast ronomical Software Applications (C ASA) 
software package (|Cornwellll2008l ; iMcMullin et alj|2007h . In 
Miriad's joint deconvolution scheme, all subfields are in- 
verted as one image, and the dirty Stokes Q and U maps are 
set to the same gridding as the total intensity image. The 



algorithm performs a maximum entropy method deconvolu- 
tion simultaneously for all polarisations, using morphologi- 
cal information from bright emission in the total intensity 
image to solve for the Q and U polarised Stokes images. 
Images are then restored using a gaussian with the same 
half-power beam width of the synthesised beanjf). While this 
method does produce a morphologically correct total inten- 
sity image, the MEM algorithm is not optimum for comput- 
ing total intensities, causing integrated intensity measure- 
ments to be systematically high by about 1 to 20 percent. 
For this reason, we used CASA's multi-scale clean (again 
joining mosaic subfields) to deconvolve images and gather 
intensity information. While CASA does compute polarisa- 
tion values, a known error in the data interface for ATCA 
polarisation information currently renders their output er- 
roneous. 

Polarisation images were created by combining aligned 
Q and U images pixel by pixel, calculating total polarised 
intensity, P, using P = \/Q 2 + U 2 . A noise term, a p was 
estimated from the RMS of a region of the Stokes Q image 
that did not contain source emission. Polarised signal below 
3cr p was considered undetected and was masked from the 
polarised intensity map. The masked polarised maps were 
then corrected for a debiasing factoiQ using a first order 
correction that again uses a p : P final = \JP 2 — (J%. 

Position angle maps are also calculated on a per pixel 
basis, with ip = | arctan^, where ip is the observed angle 
of polarisation. Vector maps of fractional polarised inten- 
sity were formed by dividing the masked total polarisation 
maps by the total intensity maps; in this way, the maps 
will only show fractional polarisation above the detection 
threshold. Polarisation vectors were then calculated using 
lengths determined by fractional polarisation levels, and po- 
sition angles from PA maps. Note that the orientation of a 
polarisation "vector" is ambiguous by 180 degrees, and does 
not provide a preferred electric or magnetic field direction 
along the printed rods. 



4 IMAGES AND DISCUSSION 

4.1 Results and Image Description 

Table [2] gives various observed and calculated parameters 
from the data. The columns are as follows: (1) Source name; 
sources with a detected nucleus have a core position fit- 
ted with a point source model, which also provides a mea- 
surement for the flux density of the nucleus; (2-3) J2000 
position of galaxy nucleus, if detected; (4) Core flux den- 
sity and error in Jy. (5) Total flux density and error in Jy; 
sources with lower limits for total integrated flux density 
are those for which we believe we have acquired a minimal 
measurement for the flux density due to observing a sub- 
region of the source. Total source fluxes were taken from 



2 See http://www.atnf.csiro.au/computing/software/miriad for 
further details of mosaic observations 

3 The debiasing term is necessary to correct the Ricean bias that 
occurs from the calculation of P from Q and U terms; there are 
multip le methods of correction as detailed in iLeahv and Ferninil 
dl989T> 
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Source Name 


RA (core) 


Dec (core) 


Score 


Sl8 


Pl8 


ni 8 


a« 


§23 


PKS 0131-36 


01:33:57.9±3.0 


-36:29:35.3±1.9 


0.03±0.01 


>0.44 


0.02* 


4.6 






Pictor A 


05:19:49. 7±0. 8 


-45:46:43.7±0.5 


1.32±0.04 


6.32±0.11 


0.50±0.06 


7.9 


-0.70 


6.80 


Centaurus A 


13:25:27.6±0.7 


-43:01:04.9±0.4 


5.98±0.17 


>28.35 


3.81* 


13.4 




46.2t 


PKS 1333-33 


13:36:39. 0±4.1 


-33:57:57.7±2.6 


0.30±0.05 


>0.74 


0.07* 


9.8 




1.70 


Centaurus B 


13:46:49. 1±0.3 


-60:24:30.0±0.2 


5.02±0.06 


8.89±0.43 


0.08±0.01 


0.9 


-0.87 




PKS 1610-60 


16:15:05.6±8.8 


-60:54:27.1±8.9 


0.14±0.05 


2.11±0.04 


0.128±0.008 


6.0 


-0.93 


1.7t 


PKS 2153-69 








3.40±0.21 


0.05±0.03 


1.5 


-0.96 


3.60 


PKS 2356-61 


23:59:04.9±7.8 


-60:55:03.4±6.1 


0.09±0.03 


1.64±0.05 


0.032±0.004 


1.9 


-0.94 


1.80 



Table 2. See section 14.11 for more detailed column descriptions: (1) Common source name (2-3) J2000 RA/Dec (4,5) Core and total 
integrated flux density and error in Jy (6) Integrated polarised intensity and error in Jy. * indicates the integrated polarisation over an 
observed subregion (7) Fractional polarisation (100 ■ Pis/Sig) ( 8) 4.85-18GHz spe ctral index (9) Flux density as detected in WMAP 
maps, with f mar king fluxes taken from the NEWPS catalogue l lLopez et a.1.1 [2007? ) . Note that the Centaurus A measurements can be 
taken as actual values (not simply limits) when referring to the inner lobes of Centaurus A. 



an integration over the MSclean-deconvolved maps, nor- 
malised by the synthesised beam. Note that this process 
has corrected for the primary beam. Errors were quanti- 
fied by the image RMS; (6) Integrated detected (> 3<r p ) 
polarised intensity in Jy (P = a/ (53 Qi) 2 + (53 Ui) 2 , where 
the sums represent the integrated and beam-normalized val- 
ues over the P > 3a p and Stokes I > 3<r pixels in the 
deconvolved Stokes Q and U images). Because of possible 
depolarisation effects when imaging the missed source re- 
gions, we cannot provide accurate total integrated polari- 
sation measurements for sources which were incompletely 
observed. Such measurements, indicated with an asterisk, 
give the integrated value of the observed source subregion; 
(7) The fractional polarisation of the observed source re- 
gion in percentage (8) The spectral index between 4.850 and 
18.0 GHz, using PMN flux densities for all so urces except 
Pictor A, where the 5 GHz measurement of iKuehr et al.l 
|l98ll ) was use d and Centaurus B, for wh ich the 5 GHz 
measurement of I Wright and Otrupcekl <ll990l) w as used ; (9) 
WMAP flux density at 23 GHz ijBennett et al.ll2003l ); flux 
densit ies marked with f are as given in the NEWPS cata- 
logue ijLopez et al.ll2007l) . 

The total intensity and fractional polarisation vector 
maps whose production are described in the previous section 
are presented in Figures OH Greyscale images show the to- 
tal intensity of emission with an "X" indicating the position 
of the suspected host galaxy for each source based on previ- 
ous identification (as noted in section with an associated 
optical galaxy. Sources for which we did not image the full 
extent of the low-frequency radio structure are accompanied 
by a radio map from previous low-frequency observations to 
give a visual reference of the observed subregion. Greyscale 
images have a linear colour gradient. 

The polarisation images indicate the observed degree 
and angle of polarisation E-field vectors, with contour levels 
as indicated below the image. The synthesised beam and a 
vector showing 100 percent polarisation are indicated in the 
lower corner of the image. 



4.2 General Properties at 20 GHz 

There are several striking and perhaps expected similari- 
ties between all the sources. Six of the objects (apart from 
PKS 2153-69) have clearly detected emission from their nu- 
clei, which in some cases were not easily visible in lower- 



frequency images; all cores have either undetected or very 
low values of polarisation (an upper limit on the polarisation 
of the cores can be calculated by taking 3 ■ op from Table [2}. 
The undetected polarisation of the cores could be the result 
of several different effects. The low values of polarisation can 
be explained by a complex scattering medium surrounding 
the inner regions of jet formation, or by a dense plasma. 
Furthering this supposition are the high (> 1000 rad m~ 2 ) 
rotation measures noted in the inner regions of some AGN. 
However, studies of galactic nuclei have been yet unable to 
characterise the spectral properties for emission from ac- 
tive galactic nuclei; the complex spectra are unable to be 
fit by pure synchrotron self absorption or free-free spectra. 
An equally likely case is the beam depolarisation of close, 
compact regions of structure that we have not resolved due 
to our lower angular resolution. With milliarcsecond reso- 
lution, the p il ot VL BA Imaging Polarimetry Survey (VIPS, 
iTavlor et al.l (|2005l )) has demonstrated that a sample of 24 
active galactic nuclei average below ~3 percent polarisation 
at 15 GHz, however the full VIPS sample will provide a 
more numerically significant assessment for high resolution 
polarisation. 



In several of the sources (esp. in the lobes of Pic- 
tor A, Centaurus A, and the inner hotspots of PKS 1610- 
60), a notable edge- brightening effect in fractional po- 
larisation is visible. While this is a feature common to 
many extended extragalactic sources, it appears in only 
a few of the sources here. The interp retation given by 
iPerlev. Roeser. and M eiscnh eimerl ljl997l ) for this effect and 
of the perpendicularity of polarisation alignment to isocon- 
tours of total flux density in the western hotspot in Pictor A 
calls upon the change in field geometry across the lobe. This 
interpretation can be supported by the absence of these ef- 
fects in the geometry and strength of observed polarisation 
in the linear jet regions in several of the sources here (for 
instance, the narrow regions of PKS 1610-60 at a ~ 16:15:27 
and a ~ 16:14:35-50, and of PKS 2356-61 at a ~ 23:59:17 
and a ~ 23:58:53). In these cases, if there is insignificant or 
no rotation of polarised signal, the magnetic field along the 
jets appears to be aligned with the jet axis, and across the 
entire region there is no apparent change in geometry, nor 
any apparent limb brightening in polarisation. 
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Figure 1. The spectrum of integrated intensity of Pictor A, in- 
cluding measurements from selected literature which have mea- 
sured its integrated flux density. Overplo tted is the spectral fit of 
IPerlev. Roeser. and Meisenheimerl l|l997l ). 



4.3 Individual properties at 20 GHz 

The inner jets of PKS 1333-33 are unique for this sample in 
their alignment; the perpendicular change in magnetic field 
orientation witnessed in transit ion from the core to the jets 
agrees with that observed by iKilleen. Bicknell. and Ekersl 
(1986). The reason for this transition is unclear, though 
could be explained through relativistic shock interactions 
with a medium surrounding the nucleus. Transverse orien- 
tations of magnetic fields are often observed in BL Lac- 
ertae objects, for in stance in the VLBI observations of 
iGabuzda etatl (|2000l ). 

Centaurus B in addition shows a striking alignment fi- 
delity in the outer regions of the source, with magnetic fields 
running perpendicularly to the jet axis along the brighter ar- 
eas of the lobes. There is a sharp change in direction and a 
brightening at the edge of the lobes, and the core is polarised 
at approximately 90 degrees to the jet emission. 

In a qualitative way, we can compare the observed elec- 
tric field vectors in our images with the lower frequency im- 
ages of polarised structure published in previous papers that 
have been corrected for rotation measure. Often images are 
published of magnetic field orientation, in which case vec- 
tor images will show vectors perpendicular to the electric 
field lines. In the case of Pictor A, th e 5 GHz images of 
IPerlev. Roeser. and Meisenh eimerl (|l997l ). which have been 
corrected for a rotation measure of 45 rad/m 2 (which cor- 
responds to a rotation at 18 GHz of 0.7 degrees) are in 
good alignment with our uncorrected images, where polar- 
isation has been detected by both observations across all 
re gions of the source. The deta i led 20 cm images published 
bv lKilleen. Bicknell. and Ekersl (1 19861 ) of the inner lobes of 
PKS 1333-33 show rotation-corrected magnetic fields in the 
lobes that are perpendicular to the jet axis. This is again in 
agreement with the results of our polarimetry, which show 
electric field vectors running directly parallel to the jet axis. 
To quantify the relative changes in spectral index and polari- 
sation for different source components in each source, a more 
rigourous wide field multifrequency observing campaign will 
be necessary. 



4.4 Prospects for CMB mission calibration 

Beam depolarisation imposes a restriction on sources as suit- 
able CMB mission calibrator candidates. The requirements 
for a good flux density and polarisation calibrator for the 
lower end of Planck's observing frequencies (the lowest band 
centre on the LFI is 30 GHz) include integrated emission > 
1 Jy across the observation bands, polarised intensity at lev- 
els > 200 mjy, and unresolved emission within the telescope 
beam. The LFI b eam resolutions range from ^33 ' at 30 GHz 
to 14' at 70 GHz (|The Planck Consortium! 120051 ). The flight 
path and pointing path of the Planck Satellite are such that 
the north and south ecliptic poles will be observed often, 
such that a source in the region around the south ecliptic 
pole would be in an ideal position. A limitation on source 
size coupled with the depolarisation that will certainly oc- 
cur within the large Planck beams means that the brightly 
polarised northern inner lobe of Centaurus A, for example, 
cannot be easily isolated and used for calibration. 

The best candidate to arise from this data set is Pic- 
tor A, which fits the Planck requirements in several ways. 
In total, the source extends approximately 10 arcminutes; 
the source yields ~ 0.50 Jy of integrated polarised inten- 
sity, which is dominated by the brightly polarised western 
hotspot. Because the emission from the hotspot is highly 
diffuse and many light years in linear size, total and po- 
larised emission from this region cannot change rapidly with 
time, making the source particularly useful for the long ob- 
serving term of the Planck mission and more useful than 
highly polarised sources which are more compact, variable 
sources. Conveniently, Pictor A is located within 20 degrees 
of the south ecliptic pole. The usability of Pictor A as a high- 
frequency calibrator is limited by its angular size, which may 
be marginally resolved at the highest frequency, and its spec- 
tral index, which is in decline at v > 20 GHz; both of these 
could limit the use of Pictor A as a flux density calibrator 
at the upper ends of the LFI bands. The integrated spec- 
trum of Pictor A is plot ted in figure Q] Overplotted is the 
log-lin ear spectral fit of IPerlev. Roeser. and Meisenheimerl 
(1997), which was fit to data below 5GHz; the higher fre- 
quency WMAP data and the data from our observations 
generally follow this log-linear trend, implying that this 
spectral index will continue into the relevant Planck bands. 
Currently, Pictor A stands as the best candidate from all 
known extragalactic sources. 



5 SUMMARY 

We have used interferometric wide-field imaging techniques 
to image seven of the brightest and most extended sources 
in the 20 GHz southern sky region —30° > 8, with the 
goal of enhancin g the statistical total intensity and polari- 
sation samples of iMassardi et af1 l|2008l ) and Burke et al. (in 
prep) . The data analysis included flux density and polarisa- 
tion measurements and maps for each of the seven sources. 
The images have revealed clearly detected 20 GHz emission 
from all galactic nuclei except that of PKS 2153-69. In all 
cases, the fractional polarisation of the galactic core was not 
detected or less than 1 percent (in the case of Centaurus A 
and PKS 1333-33). While edge-brightening and a perpen- 
dicularity of electric field polarisation vectors to brightness 
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isocontours was observed in some sources, the lack of such 
features in compact jet regions supported the hypothesis 
that the features may be due to reorientation of field geome- 
tries across broad lobe and hotspot components. All sources 
for which we imaged the entire source had declining spectra 
across 5 to 18 GHz at an average spectral index of -0.88. 

PKS 1333-33 was one of the three sources exhibiting de- 
tected polarisation in its core, and it showed a perpendicular 
field geometry to the close jets spanning the core. The jets 
exhibited a transverse magnetic field configuration to the 
jet axis, a feature usually observed in BL Lacertae objects. 
Centaurus B exhibited a similar change in alignment be- 
tween its core and outer regions with very little variation in 
angle along the peaks of the eastern and western hotspot re- 
gions, however the smaller-scale jets did not have detectable 
polarisation. 

We have put forward Pictor A as an excellent extra- 
galactic calibrator candidate for the ESO Planck Satel- 
lite mission that aims to measure the polarised signature 
of the CMB. Pictor A's integrated polarised intensity at 
0.50±0.06 Jy, its total intensity of 6.32±0.11 Jy, its position 
in the sky and its angular size contribute to its suitability 
as a candidate. 
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Figure 2. PKS 0131-36: The beam shape is indicated in the lower right hand corner of the images in Panels (a) and (b); (a) Radio refer- 
ence map and greyscale image of total intensity. An "X" indicates the position of the associated optical galaxy. The reference radio map was 
constructed using an 843 MHz image from the SUMSS postage stamp server, available at http://www.astrop.physics.usyd.edu.au/SUMSS, 
while the greyscale image was created from our 18 GHz mosaic data. Note that the apparent termination of the source at the edges of 
the observation field is due only to decreased sensitivity near the edge of our field of view, (b) 18 GHz polarisation map; vectors show the 
observed degree and orientation of the electric field vectors, while contours trace the total intensity of the radio source at levels given in 
the text below the image. The rod in the lower right hand corner of the image indicates the length of a 100 percent polarised component. 
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Figure 3. Pictor A: The beam shape is indicated in the lower left hand corner of the images in Panels (a) and (b); (a) The 18 GHz 
greyscale image of total intensity. An "X" indicates the position of the associated optical galaxy, (b) 18 GHz polarisation map; vectors 
show the observed degree and orientation of the electric field vectors, while contours trace the total intensity of the radio source at levels 
given in the text below the image. The rod in the lower left hand corner of the image indicates the length of a 100 percent polarised 
component. 
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Figure 4. Centaurus A: The beam shape is indicated in the lower left hand corner of the images in Panels (a) and (b); (a) Radio 
reference map and 18 GHz greyscale image of total intensity. An "X" indicates the position of the associated optical galaxy. The 5.0 GHz 
reference radio map was reproduced from Ron Ekers' 1977 single-dish observations, (b) 18 GHz polarisation map; vectors show the 
observed degree and orientation of the electric field vectors, while contours trace the total intensity of the radio source at levels given in 
the text below the image. The rod in the lower left hand corner of the image indicates the length of a 100 percent polarised component. 
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Figure 5. PKS 1333-33: The beam shape is indicated in the l ower right hand corner of the images in Panels (a) and (b). (a) 843 MHz 
radio reference map reproduced from I Jones and McAdaml (l992) with permission from Bruce McAdam, and 18 GHz greyscale image of 
total intensity. An "X" indicates the position of the associated optical galaxy, (b) 18 GHz polarisation map; vectors show the observed 
degree and orientation of the electric field vectors, while contours trace the total intensity of the radio source at levels given in the text 
below the image. The rod in the lower right hand corner of the image indicates the length of a 100 percent polarised component. 
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Figure 6. Centaurus B: The beam shape is indicated in the corner of the images in Panels (a) and (b). (a) The 843 MHz Radio 
reference map for the observed regions covering the core, eastern lobe, and western lobe, reproduced from SUMSS postage stamp server. 
(b) 18 GHz total intensity grcyscale for the eastern lobe, core, and western lobe, respectively. An "X" indicates the position of the 
associated optical galaxy. (c) IS GHz polarisation maps for the corresponding regions; vectors show the observed degree and orientation 
of the electric field vectors, while contours trace the total intensity of the radio source at levels given below the image. The rod in the 
corner of each image indicates the length of a 100 percent polarised component. 
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Figure 7. PKS 1610-60: The beam shape is indicated in the lower left hand corner of the images in Panels (a) and (b). (a) The 
18 GHz greyscale image of total intensity. An "X" indicates the position of the associated optical galaxy, (b): 18 GHz polarisation map; 
vectors show the observed degree and orientation of the electric field vectors, while contours trace the total intensity of the radio source 
at levels given in the text below the image. The rod in the lower left hand corner of the image indicates the length of a 100 percent 
polarised component. 
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Figure 8. PKS 2153-69: The beam shape is indicated in the lower left hand corner of the images in Panels (a) and (b). (a) The 
18 GHz greyscale image of total intensity. An "X" indicates the position of the associated optical galaxy, (b) 18 GHz polarisation map; 
vectors show the observed degree and orientation of the electric field vectors, while contours trace the total intensity of the radio source 
at levels given in the text below the image. The rod in the lower left hand corner of the image indicates the length of a 100 percent 
polarised component. 
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Figure 9. PKS 2356-61: The beam shape is indicated in the lower right hand corner of the images in Panels (a) and (b). (a) The 
18 GHz greyscale image of total intensity. An "X" indicates the position of the associated optical galaxy, (b) 18 GHz polarisation map; 
vectors show the observed degree and orientation of the electric field vectors, while contours trace the total intensity of the radio source 
at levels given in the text below the image. The rod in the lower right hand corner of the image indicates the length of a 100 percent 
polarised component. 
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